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ABSTRACT 

A new generation of sensitive X-ray measnrements are indicating that the 
existence of X-ray attennation colnmn densities, > 10^^cm“^ is qnite common 
amongst broad absorption line qnasars (BALQSOs). This is signihcant to the 
geometry of the broad absorption line (BAL) ontflow. In particnlar, snch an 
X-ray shield also shields eqnatorial accretion disk winds from the UV, thereby 
preventing high velocity eqnatorial ontflows from being lannched. By contrast, 
bipolar winds initiated by continnnm radiation pressnre from the fnnnel of a slim 
accretion disk flare ontward (like a trnmpet) and offer vastly different absorbing 
colnmns to the X-ray and UV emission which are emitted from distinct regions 
of the disk, ~ 6M and ~ lOM — 40M, respectively (where M is the radins of the 
black hole). Recent nnmerical work indicates that it is also possible to lannch 
bipolar ontflows from the inner regions of a thin disk. The recent discovery with 
VLBI that the Galactic analog of a BALQSO, the X-ray binary Circinns X-1 
(with high velocity P Cygni X-ray absorption lines) is viewed virtnally along the 
radio jet axis (and therefore along the spin axis of the black hole and the normal 
to the accretion disk) has rekindled interest in the bipolar models of BALQSOs. 

We explore this possibility by stndying the nearest BAL QSO, MRK 231. High 
resolntion 2-D optical spectroscopy and VLBI mappings of the radio jet axis 
indicates that the BAL ontflow is parallel to the parsec scale radio jet. 

Subject headings: qnasars: absorption lines — X-rays: galaxies — galaxies: active 
— accretion disks — black holes 
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1. Introduction 

It is of great physical significance that about 10% - 15% of quasars show broad absorption 
line (loosely dehned as absorbing gas that is blue shifted at least 5,000 km/s relative to 
the QSO rest frame and displaying a spread in velocity of at least 2,000 km/s) outflows 
(Weymann 1997). It is widely believed that all radio quiet quasars have BAL flows, but the 
designation of a quasar as a BALQSO depends on whether the line of sight intersects the 
solid angle subtended by the outflow. The standard model of quasars is one of a hot accretion 
flow onto a black hole and a surrounding torus of molecular gas (Antonucci 1993). A very 
popular model of the BALQSOs is an equatorial wind driven from the luminous disk that 
is viewed at low latitudes, just above the molecular gas (Murray et al 1995). Alternatively, 
one can have a bipolar BAL flow that conforms with the observations (Punsly 1999b; Proga 
and Kallman 2004). However, BALQSOs are so distant that direct imaging of the BAL 
region is beyond the resolution of current optical telescopes. Thus, much of the discussion of 
BAL geometry is based on deductive reasoning. In this Letter, we implement observations 
of nearby BAL systems in order to achieve the high resolution necessary to determine the 
geometry directly. In particular, we discuss new observations of the nearest BALQSO MRK 
231. 


Recently, new studies of the X-ray binary Circinus X-1 has changed our perspective. 
Circinus X-1 has P Cygni absorption prohles in the X-rays associated with a high velocity 
broad absorption line flow. Thus, the accretion disk of Circinus X-1 has been considered 
a Galactic example of a broad absorption line quasar (Brandt and Schulz 2001). Radio 
observations of Circinus X-1 have now indicated jet outflow velocities of 15 c (Fender et al 
2004). This can only occur for a relativistic flow that is viewed within a few degrees of the 
axis of the jet (Lind and Blandford 1985). The Galactic example of a broad absorption line 
quasar is therefore viewed pole-on. 

The findings on Circinus X-1 raise the question, are BAL QSOs being viewed close to 
the black hole spin axis? We begin this empirical analysis of BAL geometry by studying 
the wealth of new X-ray data on BALQSOs. A picture seems to be emerging that shows 
Thompson thick absorbing columns to the X-ray source. In the most realistic models of 
equatorial outflows, Murray et al (1995); Proga et al (2000), Thompson thick columns are 
catastrophic because they block the UV photons form accelerating the BAL region outward. 
Thus, the new X-ray observations indicate that high velocity equatorial BAL winds as we 
know them can not be attained in many BALQSOs. In order to determine the geometry of a 
true BALQSO, one must turn to the nearest example in order to achieve adequate resolution, 
MRK 231. Fortunately, MRK 231 also has a radio jet that is powerful enough to be imaged 
on pc scales. In this Letter, we combine VLBI mapping of the jet with 2-D spectroscopy of 
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MRK 231 to show that both the radio jet and the BAL region flow outward from the nucleus 
along PA ^ —120°. 


2. X-ray observations of BALQSOs 

The X-ray flux from the quasar tends to over-ionize gas the atoms in a putative BAL 
outflow, preventing resonant absorption and the resultant line driving force. In a break¬ 
through paper of Murray et al (1995), it was shown that if the quasar X-rays were shielded, 
a line driven wind was possible from the accretion disk which they argued was preferable 
to an outflow of dense BAL clouds and the associated confinement problems. The earli¬ 
est X-ray observations of BALQSOs supported the notion of X-ray shielding gas. ROSAT 
soft X-ray measurements indicated that BALQSOs were X-ray suppressed relative to other 
quasars (Green and Mathur 1996). Numerically, further evidence was found in Proga et 
al (2000) in which models naturally produced a shield of X-ray emitting gas between the 
central radiation source and the high velocity winds. In this section, we consider the most 
recent X-ray observations and their application to the details of the shielding gas. 


2.1. Observing X-ray Absorption Colnmns 

The ROSAT observations of BALQSOs indicated that X-ray absorbing columns > 
lO^^cm”^ were quite typical (Green and Mathur 1996). However, as observations at higher 
energy, higher sensitivity and higher resolution have become available estimates of these 
absorbing columns have increased dramatically. This essentially is a consequence of the fact 
that at Nh = 10^^cm“^ the absorbing column is essentially black to soft X-rays (Murray et 
al 1995). One needs higher energy X-ray measurements to learn more. Secondly, the fluxes 
from the BALQSO nucleus are so attenuated that they can be swamped by any nearby X- 
ray source. Thus, high resolution telescopes such as XMM are valuable. We illustrate these 
points with a couple of examples. 

The prototypical BALQSO is PHL5200. Deep ASGA observations indicated an absorb¬ 
ing column of Nh ~ 5 x lO^^cm”^ Mathur et al (2001). However, deeper high resolution 
imaging with XMM in Brinkmann et al (2002) indicated that the poor resolution of ASGA 
allowed X-rays from a background radio source to flood the detectors. Only a small per¬ 
centage of the total X-ray flux were attributable to PHL 5200. There were not even enough 
photons detected in ~ 46 ksec to determine the spectral properties and the absorption. One 
thing that is certain, the absorbing column is clearly significantly larger than the ASGA 
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value of Nh ~ 5 x lO^^cm ^ and therefore Thompson thick. 

Next we review the X-ray data on MRK 231, the nearest BALQSO. Using ASCA and 
ROSAT observations, Turner (1999), determined an absorbing column of ~ 2 x 10^^cm“^ 
to Nh ~ lO^^cm”^ depending on whether the absorbing material is neutral or ionized. 
The ASCA data were analyzed in Maloney and Reynolds (2000) who deduced that the 
direct observation of the nuclear X-rays was blocked by an optically thick absorber and we 
observe a strong reflected and scattered component through a neutral column density of 
Nh ~ 3 X lO^^cm”^. Similarly, from Chandra data Gallagher et al (2002) hnd a best fit 
to the data if the central source is blocked by a Compton thick absorber, Nh r'U lO^^cm ^ 
and we detect scattered X-rays through multiple lines of sight through absorbing columns 
of Nh ~ lO^^cm”^ - Nh ~ 6 x lO^^cm”^. Finally, hard X-ray BeppoSAX observations of 
MRK 231 combined with XMM deep observations have been made by Braito et al (2004). 
The high energy coverage combined with the 20 ks XMM observation provides the best 
determination of spectral parameters to date. They find that Nh ~ 2 X 10^"^cm Thus, as 
the measurements have improved, the column density estimates have increased. 

We also note the deep ASCA observation of PG0946-I-301 in Mathur et al (2000). After 
previous non-detections only yielded upper limits to Nh, this deep observation was best £t 
with a Thompson thick column of Nh ^ 10^^cm“^. The large column densities that have 
been deduced from the X-ray data are clearly only lower limits. Once the absorbing column 
is black, small amounts of scattered X-rays or X-rays leaking through a patchy absorber 
(small domains with Nh ~ 10^^cm“^ to Nh ~ lO^^cm”^) will dominate the X-ray flux and 
Nh of the bulk of the absorbing column will be greatly underestimated. Therefore, in view 
of the recent X-ray data it is now apparent that many BALQSOs are likely to have X-ray 
absorbing columns with Compton optical depths considerably larger than one. 


2.2. The Implications for Disk Wind Models 

In order for an equatorial flow to be launched from an accretion disk, the wind must 
emanate at large distances from the hole, r > 10^® cm, so that the radiation source appears 
as a flux of photons roughly parallel to the accretion disk surface. This allows momentum 
transfer along the equatorial direction, driving the equatorial BAL wind. As the wind is 
launched closer and closer to the source it looks more and more bipolar, since the photon 
flux attains a large vertical component (see Proga et al (2000) and references therein). 

A second necessary ingredient for the geometry is a region of X-ray shielding gas in 
order to insure the proper ionization state is achieved at r > 10^® cm so that lithium like 
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atoms are available for the line driving force to work efficiently (Murray et al 1995). The 
discovery of the large X-ray absorbing columns in BALQSOs has been heralded as the likely 
verihcation of the shielding gas at the inner edge of the wind and a triumph of the model 
Gallagher et al (2002). However, this is not really shielding gas, but hltering gas: it must 
filter out the X-rays and allow the UV emission to pass through with small attenuation. By 
geometric considerations, the UV and X-ray emitting regions are well inside of the distant 
wind zone and the shield (Proga et al 2000). Thus, if the X-rays are shielded, so are the 
UV photons. For this reason Proga et al (2000) states explicitly that for Nfj > 10^"^cm“^, 
’’the line force form the central engine could be so much reduced that it may not be able 
to accelerate the gas to high radial velocities and no strong wind will be produced.” Thus, 
the large column densities that have been found are problematic for equatorial models of 
BALQSO winds. However, there are disk wind models of bipolar outflows that do not suffer 
these same problems. Bipolar outflows are ejected from the funnels of slim disks in Punsly 
(1999b) and the work of Proga and Kallman (2004) strongly suggest that bipolar winds can 
be driven from a thin accretion disk at the relatively small radii. 


3. The BAL Outflow Geometry of MRK 231 

Bipolar models of BALQSOs do not suffer the UV attenuation problem that arises from 
the large X-ray absorbing columns (Punsly 1999a,b). The wind is launched by continuum 
radiation pressure and subsequently boosted by line driving forces farther out. Thus, near 
its base the outflow is a hypersonic under-expanded jet that flares outward like a trumpet. 
Thus, regions at different cylindrical radii on the surface of the disk (i.e., such as the X-ray 
region, r ~ 6M — 7M, and the UV region, r ~ lOM — 40M, where M is the radius of 
a rapidly rotating central black hole) are exposed to different absorbing columns of wind 
outflow. The X-rays typically pass through Nh lO^^cm ^ to Nh ~ lO^^cm ^ and the UV 
passes through Nh ~ lO^^cm”^ (Punsly 1999b). Combining this with the observation of a 
polar BAL outflow in Circinus X-1, it is of interest to see if any BALQSO can be shown to 
have a polar BAL wind. The best possibility for success is to choose the nearest BALQSO 
in order to get adequate resolution, MRK 231. 


3.1. VLBI Observations of MRK 231 

High resolution VLBI can be used to determine a radio jet axis at Ipc and give an 
estimate of the spin axis of the black hole and the therefore the normal to the accretion disk. 
MRK 231 has been observed with VLBI by Taylor et al (1994); Ulvestad et al (1999a,b); 
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Lonsdale et al (2003). The highest resolution measurements have been made at 15 GHz in 
Ulvestad et al (1999a,b) with a resolution of about 0.9 pc. A jet like structure propagates 
at PA = —112° in one epoch and PA = —115° in another epoch. The large variability of the 
northeast component (a factor of 2.5 in 1.8 years) allows its identihcation as the radio core. 
This identification of the core is consistent with our 2-D spectroscopy. 


3.2. 2-D Optical Spectroscopy of MRK 231 

Optical spectroscopy was performed on MRK 231 with the 4.2 m William Herschel 
Telescope in Lipari et al (2005). The detector was an array of hber optics and a spectrum 
was reduced from each of the hber optic responses. Each hber optic captured light from a 
sector of the sky 0.45” in radius. The result is a 2-D mosaic of spectra that covers the central 
region of MRK 231, 10 kpc x 13 kpc in extent. 

The results of the spectroscopy that are relevant are plotted in hgure 1. We are interested 
in features representing an outhow and the blue wing of the broad Hct line provides some 
clues. The upper spectrum is from the nuclear hber optic and encompasses the region 
r < 366 pc, Hq = 75 km/s/Mpc. The lower plot is from the adjacent hber optic at PA = 
—120°, 490pc < r < 1220pc from the nucleus. Note the broad emission line feature that is 
superimposed on the blue wing of the Ha line. This feature exists in the nuclear hber optic 
and is more pronounced in the adjacent hber optic at PA = —120°. There is no evidence 
of this blue bump in the spectra of any other of the hber optics that are adjacent to the 
nucleus (the spectra at all of the hber optic positions can be found in (Lipari et al 2005)). 
The conclusion is that this is high velocity Ha emitting gas that is howing at PA Ri —120°. 
The gas is actually radiating more blue shifted photons from within r < 366 pc, but this is 
swamped by the elevated background emission. 

One can estimate the peak velocity of the blue shifted Ha by extrapolating the back¬ 
ground broad line wing. This is indicated by the dashed lines in hgure 1. The peaks of the 
emission relative to the background wing are indicated by daggers. The extrapolation of the 
continuum is somewhat arbitrary and we have chosen are best ht so that the two spectra 
have similar background blue wings. With this ht, we use the peak of Ha at A = 6839A 
from Lipari et al (1994) and the cosmological redshift of z=0.042 to hnd the relativistic 
outhow velocities associated with the blue shifted peaks of 4,729 km/s in the nucleus and 
4,506 km/sec at PA = —120° (Lightman et al 1975). Note that there is some uncertainty 
in the estimate of the PA = —120° peak due to the emergence of narrow line emission that 
contaminates the broad prohle. We ht the background blue wings in a variety of ways that 
modify the velocity estimate associated with the peak of the nuclear outhow (varies within 
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the range of 4,500 km/sec - 5,000 km/sec depending on the background £t) and for the PA 
= —120° spectra (varies from 4,300 km/sec - 4,700 km/sec, depending on the background 
fit). 

As mentioned above the narrow line emission at PA = —120° adds some uncertainty to 
the estimate of the velocity of the blue bump emitting gas. It could be responsible for the 
slight disagreement with the nuclear blue bump velocity, since the narrow lines contribute 
flux to the red side of the blue bump. If one goes to the next fiber optic out along PA 
= —120° located between 1340 pc and 2030 pc from the nucleus, the narrow line emission 
becomes very strong (Lipari et al 2005). None of the other fiber optics show this strong 
emission. The generation of a strong narrow line region is expected from jet propagation in 
the ISM (Bicknell et al 1998). In fact, a similar spectral effect has been attributed to bow 
shocks in the Circinus galaxy (Veilleux and Bland-Hawthorn 1997). Thus, it seems clear 
that the outflow, initially traced by the inner VLBI jet, extends > 1 kpc at PA —120°. 
This direction validates the variable VLBA component as the nucleus (as opposed to PA 
60°). 


3.3. The Relationship to the BAL Gas 

MRK 231 is notable for the existence of many strong low ionization broad absorption 
lines. It is particularly noteworthy that the primary Na I D doublet. He I A3889 I, Ca II K 
A3934 and Ca II H A3968 broad absorption lines occur with velocities 4,590 km/sec - 4,670 
km/sec (Lipari et al 1994; Martel 1997). There is Ho emitting gas in the nuclear region 
that also flows at 4,690 km/sec. Either this is a coincidence or the two flows are one and 
the same. The latter seems more plausible to us and therefore, we conclude that the low 
ionization BAL gas is flowing outward at PA si —120°, parallel to the projection of the jet 
axis on the sky plane to within 10°. Detailed studies of other low ionization BALQSO place 
the BAL region ~ 200 pc -700 pc from the nucleus (De Kool et al 2001, 2002). However, since 
the blue bump is the most luminous in the nuclear fiber optic, the BAL flow can be located 
anywhere within the nuclear region, r < 366 pc and need not be at such large distances. 


4. Conclusion 

In this Letter, we have discussed how the large observed X-ray absorbing columns in 
BALQSOs pose an apparent catastrophic problem for the currently realistic existing equato¬ 
rial BAL wind models. We also point out that this is not an issue for bipolar flows. Evidence 



is presented from Circinus X-1 and MRK 231 that bipolar BAL flows seem to exist. There 
is likely many sources of BAL flows in quasars. This is particularly true for IR QSOs and 
high redshift QSO which are likely to be young objects (Lipari et al 2005). Our analysis 
indicates that bipolar BAL flows could be common. 

Motivated by our hndings it has become apparent that efforts to create new numerical 
models that are capable of driving bipolar winds from accretion disks is essential to our 
understanding of BALQSOs. One promising model is the disk wind model of Proga and 
Kallman (2004). At relatively small radii, it appears that bipolar winds can be launched. 
At present, it is difficult to compare the anticipated observed properties of these inner-disk 
driven winds with funnel driven outflows since the published model is merely an example of 
a broad class of models that was not optimized to produce for bipolar winds. The existing 
model has some apparent observable properties due to the choice of a disk from (Shakura 
and Sunayev 1973) with the far UV emitting region very distant from the black hole at 
r < 150M that creates a large extreme UV (EUV) source (A = lOOA to A = lOOOA) in 
the inner regions of the disk. This source over-produces the flux in the EUV by 2 orders 
of magnitude compared to composite HST quasar spectra (Telfer et al 2002; Zheng et al 
1997). There is a large scatter in the quasar spectral shapes used in these composites, so it 
is possible that some individual spectral energy distributions could peak deep into the EUV. 
But, this is clearly not typical. At the present level of modeling, it is not clear if this large 
EUV excess is an expected characteristic of bipolar winds from the inner disk. It might be 
that the parameters of future models will not require such an EUV excess to occur. It seems 
that a high priority in this held is to reproduce these simulations with a different disk model, 
perhaps a larger central black hole mass about a rapidly spinning black hole that will remove 
the EUV excess and perhaps enhance the bipolar wind production from the inner disk. 


We would like to thank Kirk Korista, Timothy Kallman and Robert Antonucci for input 
that helped with the revision of this manuscript. The comments of an anonymous referee 
greatly improved the hnal presentation of this discussion. 
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Fig. 1.— The 2- D spectra from the William Herschel Telescope with the INTEGRAL hber 
system of the blue wing of the Ha emission from the nuclear hber optic (upper plot) and 
the adjacent hber optic at PA =—120deg (lower plot). The left axis is the hux density of 
the nucleus and the right axis is the hux density from the hber optic at PA = —120°. The 
nuclear hux is the from a region, r < 366 pc. The hber optic at PA = —120° is a distance, 
490pc < r < 1220pc from the nucleus. 







